The neutron radii of 208 Pb and neutron stars 
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A new relation between the neutron skin of a heavy nucleus and the radius of a neutron star is 
proposed: the larger the neutron skin of the nucleus the larger the radius of the star. Relativistic 
models that reproduce a variety of ground-state observables can not determine uniquely the neutron 
skin of a heavy nucleus. Thus, a large range of neutron skins is generated by supplementing the 
models with nonlinear couplings between isoscalar and isovector mesons. We illustrate how the 
correlation between the neutron skin and the radius of the star can be used to place important 
constraints on the equation of state and how it may help elucidate the existence of a phase transition 
in the interior of the neutron star. 



What determines the size of a neutron star? For spher- 
ical, static stars in hydrostatic equilibrium, the so-called 
Schwarszchild stars, the sole feature responsible for their 
size is the equation of state of neutron-rich matter. The 
skin of a heavy nucleus — a system 18 orders of magni- 
tude smaller and 55 orders of magnitudes lighter than a 
neutron star — is also composed of neutron-rich matter, 
although at a lower density. 

In a recent publication we studied the relation be- 
tween the neutron skin of 208 Pb and the non-uniform 
solid crust of a neutron star For models with a stiff 
equation of state it is energetically unfavorable to sep- 
arate uniform nuclear matter into regions of high and 
low densities. Thus models with a stiff equation of 
state predict low transition densities from non-uniform 
to uniform neutron-rich matter and consequently thinner 
crusts. The thickness of the neutron skin in 208 Pb also 
depends on the equation of state of neutron-rich matter. 
The stiffer the equation of state the thicker the neutron 
skin. Thus, an inverse relationship was established: the 
thicker the neutron skin of a heavy nucleus the lower 
the transition from non-uniform to uniform neutron-rich 
matter. 

In this work we study the relation between the neutron 
skin of a heavy nucleus and the radius of a neutron star. 
Indeed, we will show that models with thicker neutron 
skins produce neutron stars with larger radii. Such a 
study is particularly timely as it complements important 
advances in both experimental physics and observational 
astronomy. Indeed, a proposal now exists at the Jeffer- 
son Laboratory to measure the neutron radius of 208 Pb 
via parity- violating electron scattering Moreover, a 
number of improved radii-measurements on isolated neu- 
tron stars, such as Geminga f§, RX J185635-3754 ||-§, 
Vela [Hi, and CXOU 132619.7-472910.8 |[o| are now 
available. While these measurements are not yet accu- 
rate enough to set stringent limits on the equation of 
state, they represent an important first step in that di- 



rection pn . 

Our starting point will be the relativistic effective-field 
theory of Rcf. ]l2] ] supplemented with new couplings be- 
tween the isoscalar and the isovector mesons. The inter- 
acting Lagrangian density for this model is given by [EUH 
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The model contains an isodoublet nucleon field (^) in- 
teracting via the exchange of two isoscalar mesons, the 
scalar sigma (</>) and the vector omega (V M ), one isovec- 
tor meson, the rho (b M ), and the photon (A^). In addi- 
tion to meson-nucleon interactions the Lagrangian den- 
sity includes scalar and vector self-interactions. (Note 
that while the original model allows for p-meson self- 
interactions [fl2f , their phenomenological impact has been 
documented to be small so they will not be consid- 
ered in this contribution). The scalar self-interaction 
is responsible for reducing the compression modulus of 
nuclear matter from the unrealistically large value of 
K = 545 MeV HQ, all the way down to about 
K = 230 MeV. This latter value appears to be consis- 
tent with the isoscalar giant-monopole resonance (GMR) 
in 208 Pb [ p^|Jl6| . Omega-meson self-interactions have 
proven essential for the softening of the equation of state 
at high density. Indeed, without them large limiting 
masses for neutron stars (of about 2.8M Q ) are predicted, 
even for the softer models that provide a good descrip- 
tion of the giant-monopole resonance |r|. This is be- 
cause the GMR constrains the equation of state around 
saturation density but leaves the high-density behavior 
practically undetermined. Models that include omega- 
meson self-interactions soften the high-density equation 
of state to such a degree that limiting masses of 1.8Af Q 
become possible (13] . Finally, the nonlinear couplings A s 
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and A v are included to modify the density-dependence of 
the symmetry energy jl]] . 

We compute the neutron radius of 208 Pb and the ra- 
dius of a "canonical" 1.4 solar-mass neutron star for all 
the parameter sets listed in Ref. pj. One of these is 
the very successful NL3-model of Lalazissis, Konig, and 
Ring. [jlTj. The other models (S271 and Z271) were intro- 
duced in Ref. Q . All the sets have been constrained to re- 
produce three important properties of symmetric nuclear 
matter at saturation: the saturation density (1.3 fm~ ), 
the binding-energy per nucleon (—16.25 MeV), and the 
compression modulus (271 MeV). The value of the ef- 
fective nucleon mass at saturation, which is not accu- 
rately known, as well as the strength of the omega-meson 
self-coupling coupling (£) differ in the various models . 
This in turn permits a modification of the high-density 
component of the equation of state in the different mod- 
els. 

The energy density of symmetric nuclear matter can 
be computed in a mean-field approximation by solving 
the classical equations of motion for the meson fields. In 
the mean-field limit it is given by [Q 
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Note that the following definitions have been introduced: 
<J>o = ,9s0o an d Wo = g v Vo. The equation of state for 
symmetric nuclear matter is displayed on the left panel 
of Fig |l|. As advertised £, and to a lesser extent M*, 
are responsible for a softening of the equation of state at 
high density. 
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FIG. 1. Binding-energy per nucleon in symmetric nuclear 
matter (left panel) and symmetry energy (right panel) in the 
various models discussed in the text. 

However, in order to compute the equation of state for 
neutron-rich matter one must supplement the equation 



of state for symmetric nuclear matter with the symmetry 
energy. The symmetry energy, a positive-definite quan- 
tity, is imposed as a penalty on the system for upsetting 
the N = Z balance. It is given by 
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where E* F = y/kp + M* 2 and the "effective" rho-meson 
mass has been defined as 
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In this manner, the equation of state of neutron-rich 
matter may be written as 
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where the neutron excess has been defined as 
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The symmetry energy is given as a sum of two contri- 
butions. The first term in Eq. (||) represents the increase 
in the kinetic energy of the system due to the displace- 
ment of the Fermi levels of the two species (neutrons 
and protons). This contribution has been fixed by the 
properties of symmetric nuclear matter as it only de- 
pends on the nucleon effective mass M*. By itself, it 
leads to an unrealistically low value for the symmetry 
energy; for example, at saturation density this contribu- 
tion yields ~ 15 MeV, rather than the most realistic value 
of ~37 MeV. The second contribution is due to the cou- 
pling of the rho meson to an isovector- vector current that 
no longer vanishes in the N^Z system. It is by adjust- 
ing the strength of the NNp coupling constant that one 
can now fit the empirical value of the symmetry energy 
at saturation density. However, the symmetry energy at 
saturation is not well constrained experimentally. Yet 
an average of the symmetry energy at saturation density 
and the surface symmetry energy is constrained by the 
binding energy of nuclei. Thus, the following prescrip- 
tion is adopted: the value of the NNp coupling constant 
is adjusted so that all parameter sets have a symmetry 
energy of 25.7 MeV at k F = 1.15 fin" 1 Q. Following this 
prescription the symmetry energy at saturation density 
is predicted to be 37.3, 36.6, and 36.3 MeV in the NL3, 
S271, and Z271 models, respectively (for A s = A v = 0). 

The simplicity of the symmetry energy is remarkable 
indeed. The contribution from the nucleon kinetic energy 
displays a weak model dependence through the effective 
nucleon mass and this dependence disappears in the high- 
density limit. The second term in Eq. (||) is also weakly 
model dependent, at least in the A s = A v = limit. The 
reason is simple: models constrained to reproduce the 
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symmetry energy of nuclear matter at some average den- 
sity, while maintaining the effective nucleon mass within 
the "acceptable" range of 0.6 < M*/M < 0.8, yield val- 
ues for the NNp coupling constant within 15% of each 
other. The weak model dependence of the symmetry en- 
ergy can be observed in the right panel of Fig. |l]. Note 
that the high-density behavior of the symmetry energy 
is given by 
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As they stand now, the models lack enough leverage 
to significantly modify the symmetry energy. In order to 
remedy this deficiency one must rely on the two nonlin- 
ear couplings between the isoscalar and isovector mesons 
(A s and A v ). These couplings change the density de- 
pendence of the symmetry energy by modifying the rho- 
meson mass as is indicated in Eq. For example, for 
A v 7^ the high-density behavior of the omega-meson 
field becomes 
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In either case (C = or £ ^ 0) the modifications are sig- 
nificant enough to change the qualitative behavior of the 
symmetry energy; the symmetry energy now grows lin- 
early with kp oc p 1 / 3 rather than as k F . This change in 
the qualitative behavior of the symmetry energy can be 
seen in the right panel of Fig. |l} Note that with A s ^ 
or A v ^0 an adjustment of the NNp coupling constant is 
necessary to maintain the symmetry energy unchanged 
from its fixed value of 25.7 MeV at kp = 1.15 fin -1 . Fur- 
ther, the inclusion of these nonlinear terms does not affect 
the properties of symmetric nuclear matter as b M = in 
the N = Z limit. 

The procedure described above is robust in another 
important way. While our goal is to induce changes in 
the neutron radius of 208 Pb through a modification of 
the symmetry energy, we want to do so without sacrific- 
ing the success of the models in describing the binding 
energy and charge radius of 208 Pb, both of them well 
known experimentally (B.E. — 7.868 MeV and i? c h = 
5.51 fm) p8]-|20||. That this is possible may be seen in 
Fig. ^| In this figure the neutron and proton ground- 
state densities have been computed in the Z271 model 
for three different values of the nonlinear uj-p coupling 
A v . While the softening of the symmetry energy has re- 
duced the neutron skin of 208 Pb appreciably, the charge 
radius has changed by less than 0.005 fm. 



0.10 
0.09 
0.08 
0.07 
^ 0.06 

'a 

SB 0.05 - 
cl 0.04 - 
0.03 
0.02 
0.01 
0.00 



Pb 




A =0.00 [R n -R p =0.244 fm] 

a[=0.05 [R°-r[=0.211 fm] 

- A^=0.10[R°-Rp=0.178fm] 



I i i i I 
4 5 
r(fm) 



FIG. 2. Neutron and proton densities in 208 Pb in the Z271v 
model for three different values of the nonlinear uj-p coupling 
A v . In all cases the root-mean-square charge radius is pre- 
dicted to be 7? c h = 5.51 fm. 

In Fig. H we show the equivalent plot but for an ob- 
ject 55 orders of magnitude heavier than 208 Pb: a 1.4 
solar-mass neutron star. The density profile of such a 
neutron star correlates nicely with the neutron skin of 
208 Pb. Models with a softer symmetry energy tolerate re- 
gions of large central densities thereby generating smaller 
radii. 
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FIG. 3. Density profile for a M = 1.4M@ neutron star in 
the Z271v model for three different values of the nonlinear 
uj-p coupling A v . 

Finally, the radius R of a 1.4 solar-mass neutron star 
as a function of the neutron skin R n — R p in 208 Pb is 
displayed in Fig. |J for the various models described in 
the text. All neutron-star radii were computed using the 
Oppenheimer-Volkoff equations for neutron-rich matter 
in beta equilibrium. Since this figure was generated us- 
ing an equation of state for uniform matter, it may con- 
tain small errors due to an inappropriate treatment of 
the surface region. Note that whereas the Z271 model 
has been extended to include both A s (Z271s) and 
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A v ^0 (Z271v), the other two models have the a-p cou- 
pling fixed at A s = 0. The strong correlation between 
the neutron-star radius and the neutron skin in 208 Pb 
is evident: for a given parameter set R increases with 
Rn — Rp. However, as one modifies the parameter set to 
increase M* or the w-meson self-coupling £, the equa- 
tion of state becomes softer so the pressure decreases at 
high density. As a result, the radius of the star becomes 
smaller for fixed R n — R p . For example, for a neutron 
skin of R n — R p = 0.18 fm, the radius of the star varies 
from R ~ 13 km in the NL3 model all the way down to 
i?~ 11 km in the Z271v model. Thus, we conclude that 
the radius of a IAM@ neutron star is not uniquely con- 
strained by a measurement of the neutron-skin thickness 
because R n —R p depends only on the equation of state at 
or below saturation density while R is mostly sensitive to 
the equation of state at higher densities. Yet one may be 
able to combine separate measurements of R n —R p and R 
to obtain considerable information about the equation of 
state at low and high densities. For example, if R n — R p 
is relatively large while R is small, this could indicate a 
phase transition. A large R n — R p implies that the low- 
density equation of state is stiff while a small R suggests 
a soft high-density equation of state. The change from 
stiff to soft could be accompanied by a phase transition. 
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FIG. 4. Radius of a M — IAMq neutron star as a func- 
tion of the neutron-minus-proton radius in 208 Pb for the four 
parameter sets described in the text. 

In conclusion, relativistic effective field theories that 
reproduce a variety of ground-state observables have been 
used to correlate the radius of a 1.4 solar-mass neutron 
star to the neutron skin of 208 Pb. Nonlinear couplings 
between isoscalar and isovector mesons have been intro- 
duced to modify the density dependence of the symme- 
try energy. Models with a softer symmetry energy tol- 
erate larger central densities and produce systems with 
smaller radii. Thus an important correlation is revealed: 
the smaller the skin-thickness of 208 Pb the smaller the 
size of the neutron star. Yet the radius of the neutron 
star is not uniquely constrained by a measurement of 



the neutron skin in 208 Pb. This is because the 208 Pb 
measurement constraints the equation of state only for 
densities between the transition density to non-uniform 
matter and saturation density. In contrast, the radius of 
a IAMq neutron star is mostly sensitive to the equation 
of state at high density. Yet together they provide con- 
siderable information on the equation of state. If these 
combined measurements reveal a large value of the neu- 
tron skin together with a small value of the star radius, 
this may provide strong evidence in support of a phase 
transition in the interior of the neutron star. 
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